Effect of anisotropic hopping on the Bose Hubbard model phase diagram: 
strong-couphng perturbation theory on a square lattice 
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There has been a recent resurgence of experimental efforts to quantitatively determine the phase 
diagram of the Bose Hubbard model by carefully analyzing experiments with ultracold bosonic atoms 
on an optical lattice. In many realizations of these experiments, the hopping amplitudes are not 
homogeneous throughout the lattice, but instead, the lattice has an anisotropy where hopping along 
one direction is not exactly equal to hopping along a perpendicular direction. In this contribution, 
we examine how an anisotropy in the hopping matrix elements affects the Mott lobes of the Bose 
Hubbard model. For weak anisotropy, we find the phase diagram is only slightly modified when 
expressed in terms of the average hopping, while for strong anisotropy, one expects to ultimately 
see dimensional crossover effects. 

PACS numbers: 03.75.-b 03.75Lm 37.10.Jk 67.85.Hj 
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I. INTRODUCTION 



Recent efforts in ultracold atomic gases have been 
looking at the prospects for building an optical lattice 
emulator — namely an analog quantum simulator that can 
model simplified many-body Hamiltonians used in con- 
densed matter physics, and extract important physical 
properties of the models, like their phase diagrams MQ- 
In particular, the recent experiments of Spielman [3[, et 
at, determined the critical interaction strength for the 
loss of superfluidity in the Bose Hubbard model on a 
square lattice. Their measured value is within 6% of 
the theoretically calculated values for the critical point 
{Uc ~ 15.8(20) in the experiment at a finite temperature 
versus Uc « 16.74 in theoretical calculations at T = 0). 

The optical lattice formed in this experiment was con- 
structed by a single laser beam that was reflected to 
form standing waves in two different directions, creating 
a square lattice. Due to geometrical constraints in the ex- 
perimental set up, the hopping matrix elements differed 
in the x and y directions (due to slightly different po- 
tential well depths in the two directions). The question 
we want to answer here is how much does the location 
of the critical point depend on the hopping anisotropy? 
Obviously, if the anisotropy is so large that the system 
becomes weakly coupled one-dimensional chains, then we 
expect to see a dimensional crossover of the critical be- 
havior at the tip of the lobe; in this case it would change 
from the power law behavior in two dimensions to a 
Kosterlitz-Thouless behavior in one dimension. Having a 
solid quantitative understanding of the importance of the 
anisotropy is critical for analyzing any given experiment. 

Since the critical U values are large for the Mott insu- 
lator to superfluid transition, a strong-coupling analysis 
is the most appropriate way to proceed. Strong cou- 
pling expansions have been shown to be very accurate 
in determining the phase diagram 0, i, i, especially 
when they are coupled with a scaling analysis of the tip 



of the Mott phase lobe. It is a straightforward exercise 
to modify the strong-coupling formulas to take into ac- 
count the anisotropy of the hopping matrix elements and 
determine their effect on the phase diagram. In this con- 
tribution, we will examine the effect of anisotropy on the 
two-dimensional square lattice. This is the first nontriv- 
ial system where anisotropy plays a role, and hence it 
should be the simplest to analyze as an explicit example. 
The other reason to look at this case is that experiments 
have already been performed, so we can quantitatively 
determine what the effects of anisotropy are on those ex- 
periments. 

From a heuristic standpoint, it is obvious that we 
should express results in terms of the average hopping 
matrix element i = {t^ + ty)/2, because only the aver- 
age matrix element enters into the first-order shift of the 
energies of the particle or hole defect phases of the Mott 
insulator (see below). Hence the deviations due to an 
anisotropy, which will depend on At — (t^ — ty) enter 
first quadratically in (Ai)^. This result follows from the 
fact that we should not be able to tell the difference of 
a system with a slightly larger hopping in the x direc- 
tion versus one with a slightly smaller hopping in the 
X direction, because the two cases can be mapped onto 
each other by rotating the lattice by 90 degrees (if we al- 
ways use the hopping in the y direction as our reference 
hopping). So an anisotropy of 15% should produce cor- 
rections on the order of 2%. Such corrections are likely 
to be smaller than other experimental errors, so the ef- 
fect of anisotropy is expected to be minor for realistic 
experimental situations. Indeed, this is the main conclu- 
sion of this work, but we actually find that the tip of 
the Mott phase lobe is even less sensitive to the effects 
of anisotropy so the correction factor is much smaller. 

In Section II, we derive the expressions for the strong- 
coupling analysis of the Mott insulator and of the particle 
and hole defect phases and show how the phase bound- 
aries are modified by the anisotropy. In Section III, we 
present numerical results and a scaling analysis of the 
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behavior of the critical point and how it changes with 
anisotropy. Conclusions follow in Section IV. 



II. FORMALISM 



We begin with the Bose Hubbard model whose Hamil- 
tonian is Q 

n = -Y.^^Ab, - mE"- + - 1), (1) 



where bi is the boson annihilation operator at site i, 
n.i = h\h^ is the bosonic number operator at site i, — is 
the hopping matrix element between the site i and j , U is 
the strength of the on-site repulsion, and fx is the chem- 
ical potential. The hopping matrix is nonzero for near- 
est neighbors on a two-dimensional square lattice. The 
hopping is allowed to be anisotropic, so that — ta-^^x 
corresponding to hopping parallel to the x-axis need not 
be equal to ty — tu+y, corresponding to hopping paral- 
lel to the y-axis. The lattice does maintain translation 
invariance though, so all hopping matrix elements along 
the x-direction are equal to tx and similarly for the y- 
direction. It is more convenient to express the results in 
terms of the average hopping matrix element 



t = 



H~ 



and the relative hopping matrix element 



— tx ty- 



Using these definitions it is easy to see that 



and 



(2) 



(3) 



(4) 



(5) 



These relations will become important below. 

We will be examining the stability of the Mott insu- 
lating phase with no bosons per lattice site. Following 



Ref. we need to evaluate the perturbative energy for 
adding one particle or for adding one hole into the Mott 
phase. To zeroth order in the hopping, the ground state is 
A^-fold degenerate, with N the number of lattice sites on 
the lattice. The hopping matrix breaks this degeneracy 
(to first order in the hopping). We need to find the low- 
est eigenvalue and eigenvector of the single particle ma- 
trix Sij = —tij to perform the degenerate perturbation 
theory. Since the lattice is translationally invariant, the 
lowest eigenvector is the vector (1, 1, 1, l)/\//V with 
an eigenvalue Amin = ^"^tx ^ '^ty ~ — 4i. The chemical 
potential within the Mott phase is written as the devia- 
tion from integral multiples of U: I-J, — {hq + 5)U . The 
point fi = noU corresponds to the upper tip of the hq 
Mott phase lobe, which meets the lower tip of the no -I- 1 
Mott phase lobe. The variable 6 is negative for a particle 
defect and positive for a hole defect. 

Under the assumption that the compressibility van- 
ishes continuously at the Mott phase lobe boundary, the 
transition between the incompressible (Mott) phase and 
the compressible (superfluid) phase is determined when 
the energy difference between the Mott insulator and the 
defect state (particle or hole) vanishes. The two branches 
of the Mott-phase boundary meet when 



5(P"*^(no) + l = 5(''°'°)(no). 



(6) 



The additional one on the left hand side arises because 
6 is measured from the point fi/U = no, and the hole 
defect for the no Mott insulator starts at n/U = no — 1. 
Equation ^ is used to estimate the critical value of the 
hopping matrix element beyond which the Mott-insulator 
phase ceases to exist. 

The energy of the Mott phase satisfies [1, 
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4P + (A^)^ 
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no(no + 1), 



(7) 



through third order in the hopping matrix element. No- 
tice that we used the definition in Eq. ^ and that no odd 
powers in the hopping enter because the square lattice is 
bipartite. 

Similarlv, the energy to add a particle to the system 
satisfies 
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64(2no + l)r + (25no + U)t[At^ + {Aty] - 2(2no + l)[Af^ + 3i(At)^] 
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while the energy to add a hole is 
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(9) 



Both results hold through third-order perturbation the- 
ory in the hopping. We take each equation, find the point 
where the defect energy is equal to the Mott energy (i.e., 
the LHS vanishes), and then solve for the critical chemi- 
cal potential 6 as a function of i and At. 

Note that we can obtain an estimate of the value of 
At/t needed for the dimensional crossover because the 
Mott phase lobes in one dimension have an opposite sign 
curvature from that in two dimensions, due to the cusp- 
like behavior of the Kosterlitz-Thouless transition. So a 
simple estimate for the dimensional crossover is when the 
coefficient of the quadratic term in the hopping vanishes 
for the energy differences of the defect phases. When 
no = 1, this occurs at ty/tx = 8 — V63 ~ 0.0627 for 
the particle branch of the lobe and ty/t^ = 2 — VS w 
0.2679 for the hole branch of the lobe, while the change 



in curvature occurs at ty/tx =4 — vl5 ~ 0.1270 for both 
branches when we take the limit uq —i- oo. This implies 
that we can crudely estimate the dimensional crossover 
to occur when ty/tx is on the order of 0.1-0.05, or At/t « 
1.6-1.8. 

We will use scaling theory to extrapolate the short- 
term power series expansion into a functional form ap- 
propriate for the Mott phase lobe. To begin, we define 
the parameter x by 



2t 

c7' 



(10) 



and then express the boundaries of the Mott phase lobe 
branches via a power series that includes the known 
power-law critical behavior of the tip of the lobe Q 
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A + Bx + Cx'^ + Dx^ + ... 
ix^-xf''{a + f3x + -fx^ + ...), (11) 
A + Bx + Cx^ + Dx^ + ... 



- {xc - xf" {a + /3x + -fx'^ + ...), (12) 

where the constants A, B, C, D, ... and a, /?, 7, ... can 
depend on At/i and Zv is the critical exponent which 
determines the shape of the Mott phase lobe near the 
critical point Xc. We immediately find 



A = 

B = 

C = 

D = 



1 
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(13) 
(14) 

(15) 
(16) 



The expression for the critical behavior is somewhat more 
complicated; 

{xc - xf^ia + Px + jx^) = i - x{272o + 1) 
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Note that the equality in Eq. (fTT)) is meant only in the 
sense that it holds term by term in when the LHS is 
expanded in a power series in x for small x. 

Performing the Taylor series expansion of the power 
law for small x and equating terms (keeping the series up 
to only the 7 term) produces the following four equations 
that the coefficients, Xc and Zv satisfy: 



1 (1 

2 [a 



Zvl3 



Zv 

X,. a 



Zv 



Zv{Zv - 1 



aXc 2x'^ 
Zv-i Zv{Zv-l)(i Zv{Zv -l){Zv -2) 



2ax1 



6xi? 



-(2no + l); 



6nr, -f 6no — 1 — 5nr. + 5no + 



1 



no (no + l)(2no + 1) 
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At 
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(18) 
(19) 

(20) 
(21) 



There are two ways to perform the scaling analysis. In an unconstrained analysis, we set 7 = and solve Eqs. (|18l - 
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I2ip for Zv^ Xc, cx, and /3. In a constrained analysis, we 
assume the three-dimensional XY model critical behavior 
holds for all At/i values of interest, and so we set Ziy = 
0.67 and fit Xc, a, /3, and 7. Both of these techniques 
will be compared numerically in the next section. 



III. RESULTS 



We focus on the uq — 1 Mott phase lobe since this has 
been the most accurately measured phase diagram of the 
model. For the unconstrained scaling analysis, we set 
7 = and solve for Zv, Xc, a, and (3. Using a = l/2a;f'^, 
/3 = [—6 + {Zi'/xc)]a, and defining y — Zv/xc produces 
the following two equations to sclf-consistently solve for 
the critical point and exponent 



y 



1 

Xr 



6 



1 



W(12-Xe)2+44- 42 



(22) 



(23) 
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These equations are most simply solved by iteration. 
Start with Xc — 0.1 and iteratively solve Eq. (221) for y 
and then Eq. for Xc until they stop changing. A plot 
of the relative change of the critical point Xc{At/t)/xc{0) 
versus the anisotropy At/t is shown in Fig. [1] with the 
solid line. We find the critical point initially increases un- 
til At « 0.6i where it begins to decrease. But, the effect 
is much smaller than one would have naively guessed. For 
example, if A</< = 0.1, we find Xc{At/t)/Xc{0) = 1.0004, 
which is 0.04 times smaller than the crude estimate that 
the deviation is quadratic in the anisotropy. 

These equations can no longer be solved once At > 
1.02i, where the critical Xc is pushed to zero and the ex- 
ponent vanishes. Obviously the scaling approach cannot 
faithfully determine the Mott phase lobe when At/t be- 
comes too large. In actual experiments, we expect the 
anisotropy to be small (on the order of 15% or less) so 
this is not a big issue. Surprisingly, the critical exponent 
Zv decreases with anisotropy. One might have expected 
it to increase to approach the cusp-like behavior of the 
one-dimensional system, but it doesn't do that. Perhaps 
this is because one needs higher order terms in the se- 
ries to properly estimate the exponent [6]. Note that 
at At/t = 0, we have Xc = 0.1121 and Zv = 0.5826, 
while the true critical point and exponent satisfy @, Q 
Xc = 0.1195 and Zi^ — 0.67. While this truncated analy- 
sis is accurate to about 7% for the critical point in the ab- 
sence of anisotropy, we expect the truncated analysis to 
be much more accurate for the relative change of the crit- 
ical point with anisotropy, especially for small anisotropy. 
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FIG. 1: (Color online.) Plot of the relative change of the lo- 
cation of the tip of the no — 1 Mott phase lobe as a function 
of the anisotropy. The solid (black) line is for the uncon- 
strained scaling approach and the dashed (red) line is for the 
constrained scaling approach. Note how the two curves agree 
for small anisotropy and how the effect is quite small. When 
the anisotropy becomes large enough, the series is to small to 
carry out the scaling analysis and the results become inaccu- 
rate. 



Indeed, the increase of the critical x for small anisotropy 
is precisely what we would expect, since the Mott phase 
lobe in one dimension has a cusplike shape and is pushed 
to larger x values than in two dimensions, so we expect 
the anisotropy to move the critical point to the right in 
the phase diagram. 

Next, we examine the constrained scaling theory anal- 
ysis, where we fix Zv — 0.67 and we solve for Xc, a, /?, 
and 7. In this case, we have a cubic equation to solve 
for Xc which is simple to write down by eliminating a, (3 
and 7 from the coupled equations for the critical point. 
The solution of these equations is also plotted in Fig. [T] 
with the dashed line. Note how the critical value of Xc 
always increases when we constrain the exponent; the 
method fails for At/i > 1.75. The constrained analysis 
is much more accurate for the critical point, predicting 
it to lie at Xc = 0.1173 which is only a 2% error. The 
change of Xc with anisotropy continues to be small, but 
not as small as in the unconstrained analysis; we find 
Xc{At/i)/xc{0) = 1.0010 at At/t = 0.1, which is still a 
factor of ten smaller than we originally estimated. But, 
the trend is once again to increase the critical x value, 
which is just what we expect when we compare the one- 
dimensional and two-dimensional Mott phase lobes. 

It appears that the constrained extrapolation proce- 
dure estimates the change in the anisotropy better than 
the unconstrained approach. However, in both cases, the 
change in the critical point due to anisotropy is much 



5 



too small to be observed in experiment for anisotropics 
on the order of 15% or less. But, it is imperative that 
one uses the average hopping as the energy scale in order 
for the corrections to be small. If any other relative unit 
for hopping is chosen, the changes in the critical point 
can be very large (linear in the anisotropy). 

IV. CONCLUSIONS AND DISCUSSION 

We have examined the effect of hopping anisotropy on 
the shape of the Mott phase lobe of a two-dimensional 
square lattice. When results are expressed in terms of 
the average hopping, we find the critical point changes 
by an amount that is too small to be seen experimentally 
for reasonable experimental parameters (an anisotropy of 
15% say). This implies that one need not worry about 
controlling the lattice parameters to be identical when 
examining the Mott phase transition in the Bose Hub- 
bard model. 



One can ask what will happen in three dimensions? 
Once again, if we express results in terms of the average 
hopping, we expect the errors to scale quadratically in 
the anisotropy. Furthermore, since the shape of the Mott 
phase lobe for two dimensions and for three dimensions is 
quite similar, we anticipate that the effect of anisotropy 
would be even smaller for the three dimensional case. 
Hence, we do not believe the effect of anisotropy will lead 
to any significant modification of experimental results on 
these systems as long as the anisotropy is kept at the 
modest level of about 15% or less. 
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